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Abstract. We present the preliminary analysis of the GRB light curves obtained by 
SwiftlXKi between November 2004 and December 2010. 
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^ . 1. Introduction 

c^ . 

The Swift satellite dGehrels et al. 1120041) . launched on November 2004, during the past six years 
of operation detected and observed more than 600 Gamma-Ray Bursts (GRBs). The vast ma- 
^ ■ jority (~67%) of t hese bursts were monitored in the soft X-ray band by the X-Ray Telescope 

t^ ' (XRT. iBurrows et al. 2005) starting as early as ~80 s after the trigger The standard model ex- 

t~^ , plains the X-ray afterglow of GRBs as synchrotron radiation arising from the deceleration of a 

OO ■ relativistic blast wave into the external medium. The XRT follow up, therefore, is comprehensive 

Cn I of the tail of the prompt emission and of the afterglow. The XRT sample is now large enough 

r — . to justify a statistical approach aimed at collecting, classifying and understanding the observa- 

f^ ' tional information of a wide and homogeneous sample of GRBs. This work will provide the most 

complete view of the X-ray properties of GRBs to date any existing theoretical model is asked 
to explain, while serving as a guide for future theoretical developments. The cataloguq] is now 
under completion, here we report the status of the work. 
• 5—1 

;_| ' 2. Sample, Data Reduction and Analysis 

We analysed all the GRBs detected until the end of 2010, for which the afterglow had been 
observed by XRT with enough photons to extract a measurable spectrum. The sample consists of 
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437 GRBs out of a total 658 GRBs detected by Swift of which 165 GRBs with redshift, 414 long 
GRBs (153 with z) and 23 short GRBs (12 with z). The original XRT data have been reduced 
with the method reported in lMargutti et alj (|2010|) . We extracted the XRT light curves in the 0.3- 
10 keV XRT band as well as in the 0.3-1 keV, 1-2 keV, 2-3 keV and 3-10 keV band. Our XRT 
archive contains light curves in count-rate, flux and luminosity (for the redshift subsample) and 
all the relevant parameters of interest. The light curves in flux units are calibrated accounting for 
spectral evolution and the spectra have been derived using the NH column density estimated in a 
time interval where no spectral evolution is apparent. 

3. Data analysis 

We fitted the light curves in flux units, using four functions: 



Single power-law: 

F{t)^Nt-". 



Smoothed broken power-law: 

F(t) = N 
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Sum of power-law and smoothed broken power-law: 

F(t) = A^i r"' + N2 
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- Sum of two smoothed broken power-laws: 
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where a is the decay power-law index, tb the break time, s the smoothness parameter and A^ 
the normalization. The best fit parameters were determined using the IDL Levenberg-Marquard 
least-squares fit routine (MPFIT) supplied by Markwardt (2009). All variability and fluctuations 
superimposed on the basic underlying light curve have been subtracted by iteration following 
JMargutti et al. II2OIII) . Using the fit parameters, we calculated the total fluence (and the energy 
for the subsample of GRBs with known z) of our light curves as well as the one of different parts 
of the light curves (E\, E2, £3, £4), as shown in Fig. 1; moreover we calculate the fluence (energy 
when possible) of the excesses. 



4. Classification and morphology 

We classified the XRT light curves according to the fit funct ion used, the presence of flares 
and the completeness of the curve (Fig. 2 and Tab. 1; see also iBernardini et al. 11201 1|) . For the 
different shapes (Fig. 1), we have: 

- Type 0: no breaks, Eq.[T] 

- Type la: single break, Eq.|2]and s < 0; 

- Type lb: single break, Eq.|2]and s > 0; 

- Type Ila: double broken power-law (steep-to-flat), Eq. [3j 
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Fig. 1. Cartoon of the fit function used to classify and analyse the sample of XRT light curves. We also 
indicated some of the parameters derived and the time intervals for which we measured the emitted energy. 
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Fig. 2. Classification of the light curves in the sample, as illustrated in the Tab.[T] 0, I, II, III indicate the 
fit function used where the number refers to the numbers of observed breaks in the light curve; for instance 
Type means that the light curve is a pure power law. "C" stands for complete and "U" for incomplete; 
"F" indicate the presence of flares; on the contrary, there are not flares if in the name is present a "N". For 
example, "ICN" means that the fitting function used was a smoothed broken power law with a complete set 
of the data and no flares. 



- Type lib: double broken power-law (flat-to-steep), Eq.|3j 

- Type III: double broken power-law, Eq. |4] 

A light curve is considered complete if the XRT re-pointing time is <300 s and the final count 
rate is comparable to the background. This ensures that the possible absence of the early steep 
decay is not due to an observational bias. See Tab. [T] and Fig. 1, 2 for more details about our 
sample. 



Type 


UN 


UF 


CN 


CF 


Total 





65 


4 


28 


17 


114 


I 


29 


7 


86 


31 


153 


II 


9 


4 


84 


51 


148 


III 








12 


10 


22 


Total 


103 


15 


212 


107 


437 



Table 1. Distribution of the GRBs light curves in the different classes. For more details see Fig. 2. 
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Fig. 3. The Dainotti relation (left) and the histrogram of the total X-ray energy with errors calculated with 
the Monte Carlo method (right). In red the short GRBs. 



5. Conclusion 

The temporal and spectral properties of a homogenous sample of 437 GRBs have been extracted. 
For the first time, we have a large and homogeneous data set that allows us to perform a statistical 
study of the X-ray properties of GRB afterglow. Notably, the sample of GRBs with redshift 
information comprises 165 elements: for this sub-sample it is possible to study the intrinsic 
properties of these explosions. 

From Fig. 2 and Tab. 1, we notice that the probability to observe a Type 1 or Type II XRT light 
curve is higher if the data are complete (1 17/153 for Type 1 and 135/148 for Type 11) evidencing 
an observational bias. The majority of incomplete light curves are Type (69/1 18). In addition, 
almost all the incomplete light curves have not flares ( 103/118) and 33% of the complete light 

curves have flares (107/319; see lChincarini et aini2010l) . 

As an exam ple of the output of our analysis we plot the Dainotti relation (iDainotti et al.ll20()8l 
I2OIOII2OIII) . that is the relation between the end plateau luminosity and time, as we derive it from 
the row output and the distribution (histogram) of the total X-ray energy where the uncertainties 
for each bin have been computed by Monte Carlo simulations (Fig. 3). 
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